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Quasars in the 4D Eigenvector 1 Context: 

A STROLL DOWN MEMORY LANE 

Jack W. Sulentic^, Paola Marziani^ 


Abstract 

Recently some pessimism has been expressed about our lack of progress in understanding 
quasars over the 50+ year since their discovery. It is worthwhile to look back at some 
of the progress that has been made - but still lies under the radar - perhaps because 
few people are working on optical/UV spectroscopy in this held. Great advances in 
understanding quasar phenomenology have emerged using eigenvector techniques. The 
4D eigenvector 1 context provides a surrogate H-R Diagram for quasars with a source main 
sequence driven by Eddington ratio convolved with line-of-sight orientation. Appreciating 
the striking differences between quasars at opposite ends of the main sequence (so-called 
population A and B sources) opens the door towards a unihed model of quasar physics, 
geometry and kinematics. We present a review of some of the progress that has been 
made over the past 15 years, and point out unsolved issues. 

Keywords: quasars: general - quasars: emission lines -quasars: supermassive black 
holes - galaxies: active 


1. Introduction 


The hrst hfty years of research on quasars has lead us to a paradigm involving 


three main (unresolved) components 1) supermassive black hole (SMBH; Zel’Dovich and 


Novikov 1965 Salpeter 1964, 2) an emitting accretion disk (AD; Shields 


and Sargent 1982) and 3) an obscuring torus (Antonucci and Miller| 1985). Component 


1978 Malkan 


3 arises from attempts to unify Type 1 and 2 sources on the basis of orientation. A 4th 
component would involve narrow line emission and sometimes radio-jets extending on 
scales of parsecs, kpc or even Mpc along the rotation axis of the disk, resolved or partially 


resolved in the nearest sources (Capetti et al 


Falcke et ah, 1998). The AD may 


contribute to broad line emission and as such may be a constituent of the broad line re- 


gion (BLR) (Chen et al. 

1989 

Dumont and Collin-Souffrin 

1990 

Eracleous and Halpern 

2003 Bon et al., 

2007 

), while Component 4 involves only narrow-line region (NLR) emis- 


sion unless the source is radio-loud. Cartoons showing the different components neatly 
aligned are likely too optimistic but we can hope it is often the case. Moving beyond this 
paradigm has not been easy leading to some expressions of frustration. Impediments to 
real progress involve lack of a clear dehnition of a quasar and the lack of any paradigm 
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in which to contextualize source commonalities and differences. For a long time then we 
have been stuck with the paradigm as definition and the assumption that all quasars (or, 
better said, active galactic nuclei AGN) are the same. We have been waiting for one of 
the self-styled ‘elite” to lead us out of the darkness. 

During the past 15 years, impatient with waiting, we have assembled a formalism 
designed to contextualize quasar diversity and identify the principal physical drivers of 
that diversity. This 4D Eigenvector 1 (4DE1) parameter space represents a surrogate 
H-R Diagram for type 1 AGN. We extend the definition of Type 1 to include sources 
that show both broad-line emission from the principal optical and UV permitted lines 
and sources that show optical Fell emission and that are therefore expected to accrete 
at a moderate to high-rate (dimensionless accretion rates > 10“^ — 10“^, Woo and Urry 


2002 Marziani et al.||2003c). 4DE1 was built upon pioneering studies of optical (Boroson 


and Green, 1992), UV (Gaskell, 1982) and X-ray (Wang et ah, 1996) spectra. By the 


year 2000 enough data and ideas were in place to introduce the 4DE1 formalism and 
the idea of two quasar populations A & B (Sulentic et ah, 2000a||b ) to emphasize the 
many differences between high and low accreting Type 1 sources. The populations may 
represent two distinct quasar classes - only one capable of radio-loudness ( |Zamhr et al 
2008 )-or opposite extremes of a single quasar “main sequence”. 


Developments in the statistical analysis before late 1999 are reviewed in Sulentic et al. 


(2000a). Here we retrace developments that came mainly afterwards, since the turn of the 


century, and that are associated with a better understanding and expansion of the original 
eigenvector 1 results. A recent work by Shen and Ho (2014) (hereafter SH14) basically 


confirms many past results and hopefully reawakens interest in the subject. The aim of 
this paper is review the main steps in the development and exploitation of 4DE1 mainly 
by the group of Sulentic and collaborators (Ipj - [^. We also review the “rediscovery” 


SH14 


and mention its followup (|Sun and Shen 

i 


2015). We maintain a loose 


paper 

chronological order ((|^ ^ [^ in an attempt to clarify which results can be considered 
established and which still lie on uncertain ground- therefore necessitating further study. 
We hnally point out some major open issues 


2. 2000: Formulation of an Eigenvector 1-based parameter space 

The 4DE1 parameters introduced in 2000 involve: 1) full width half maximum (FWHM) 
of broad H/? (H/3bc); 2) equivalent width (EW) ratio of the optical FellA4570 blue blend 
and broad H/3 (Rpeii)- The choice of equivalent width was motivated by its widespread 
availability in low z (<0.7) spectra. In more recent time we have considered Rpeii as 
defined from the intensity or flux ratio in order to avoid division by a continuum that is 
often steeply rising toward the blue; 3) profile shift at half maximum of high ionization 
ClvA1549 c(l/2) and 4) soft X-ray photon index (Tgoft)- The parameters are thought to 
be: 1) a measure of virialized motions in a low-ionization line emitting AD or flattened 
system of clouds that is considered an important virial estimator of black hole mass for 
large samples of quasars. 2) A sensitive diagnostic of ionization, and column density 
in BLR gas arising, as far as we can tell from shielded or outer parts of the BLR. The 
strength of optical Eell emission in many sources has been argued to support the AD ori 
gin for the emission (Gollin-Souffrin et ah, 1988 Dultzin-Hacyan et ah, 1999 poly et al. 
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2008) and a role for metallicity. 3) A strong diagnostic of winds/ontflows in the higher 
ionization broad line gas. 4) a diagnostic measnre of thermal emission likely connected 


with the accretion disk, and to the accretion state (e.g., Mineshige et ah, 2000, althongh 
see Done et al. 2012 for a dissenting view). These measnres were chosen becanse they 


were available for: 1) large nnmbers (100+) of low ^ qnasars with 2) high S/N spectra 
(S/N>20) and 3) becanse they showed statistically signihcant dispersion along the 4DE1 
main seqnence. They are Eigenvector 1 becanse they are strongly correlated and if pa¬ 
rameters 3 and 4 had been inclnded in the hrst PCA analysis they wonld have contribnted 
mnch of the power of Eigenvector 1. They are “orthogonal” in the sense that they involve 
parameters describing independent aspects of qnasar phenomenology as well as different 
physical processes connected to the BLR. 

One can mark 1998 as the year when theoretical attempts to model the BLR (of 


well stndied NGC5548) were acknowledged to have failed at least in part (Dnmont et al. 


1998). We think the pessimism that followed is reflected in the incompatibility of measnres 


for popnlation A and B qnasars (e.g. signihcantly different ionization parameter U and 
geometric factor /). After 1990, there has been a widening void between theory and 
observations, probably becanse the large diversity of qnasar properties has not been taken 
into acconnt as needed. Mnch more snccess is likely to arise from attempts at modelling 
Pop. A (e.g. I Zw 1) and B (NGC 5548) or individnal spectral types (as dehned by 
Snlentic et aL||2002 ^ below) separately. Then the qnestion of nniting the two solntions 
~ or the need to invoke separate qnasar populations (perhaps driven by a critical accretion 


rate, Marziani et al. 2014) can be addressed. 


3. 2001: The Physical Drivers of 4DE1 

After we proposed the 4DE1 contextualization we began to search for the physical 
drivers. We knew that, unlike stars, quasars are very unlikely to show the same spec¬ 


troscopic properties at different viewing angles (Wills and Browne, 1986). Boroson and 


Green (1992) suggested L/LEdd as a possible physical driver by exclusion, as a “best 
guess”. Not source luminosity (an Eigenvector 2 parameter in their PGA), not simply 
orientation (the extreme range in FWHM alone H/3 precludes that) nor geometry -all 
were discussed much more. They were writing in the early 1990s using a sample of 87 
quasars so their discussions are perhaps out of date from todays perspective. More then 
two measures would likely be needed to break the degeneracy between source orientation 


and physics. The following year (Marziani et ah, 2001) we explored the physical drivers 
of source occupation along the optical “main sequence” and concluded that source ori¬ 
entation {9) convolved with the ratio of quasar luminosity to black hole mass (L/Mbhc< 


Eddington ratio, as anticipated by Boroson and Green] 1992) could describe rather well 
source occupation in the optical plane of 4DE1. (Fig. reproduces such an overlay on 
our heterogeneous but high S/N sample (see also ESO Messenger for June 2001 Snlentic 


et al.||2001 ). It was rejected by the Nature editor. The following year two of us considered 
the role of black hole mass over 8dex range (Zamanov and Marziani, 2002). This was 


accompanied by an attempt to estimate L/LEdd from L and Mbh, the latter estimate from 


line width and scaling relations (Laor, 2000 Boroson, 2002 Marziani et ah, 2003c). The 


US author of this paper explained to his collaborators (from, Italy, Mexico and Slovenia) 


3 











































that these results would remain largely unacknowledged until someone from an “elite” 
institution rediscovered them. That’s how it works in the US at least. Many times over 
the past 14 years he lamented the slowness of his elite colleagues to grasp the signihcance 
of these really quite simple results. Finally it has come to pass (SH14)! And with the 
expected endorsements and encomia. 


At this point it might be useful to compare the abstract of SH14 with that of Marziani 


et al. (2001). In one sense, SH14 is “one step forward and one step back”. The new paper 


uses a much larger sample of quasars drawn from SDSS (one step forward) but, most of 
the spectra that they use and measure (Shen et ah, 2011) are unsuitable to produce an 
equal or greater description of the role for \j/ljEdd (one step back). SDSS is a gold mine 
but at the same time it is a mineheld for the uninitiated user. Automated processing 
of these spectra is indeed a walk through a mineheld. We have emphasized since 1996 
the need for, and importance of, quasar spectra with high enough S/N and resolution to 
permit reliable spectroscopic parameter measures. That is why Boroson and Green (1992) 
was such a step forward and their measures are still used 25 years later. 

Although in some places it may be treasonous to point it out, most low S/N quasar 
spectra are not suitable for useful spectroscopic measures. If this statement were untrue 
then, in the context of future quasar studies, most of the justihcation for 10m class 
telescopes - beyond observing, sure to be proposed, obscured quasars - would be negated. 


Figure^ shows a comparison of automated measures of SDSS quasars ( 

Shen et al. 

2011) 

used in the SH14 study with IRAF SPECFIT measures of 

Zamhr et al. 

(2010) 

. The 


comparison involves the brightest quasars in SDSS-DR5 and therefore generally showing 
S/N>20. Continuum S/N not to be confused with S/N computed in regions with strong 


emission lines. Surprisingly many of the FWHM H/d measures used in this paper (Shen 


et al., 2011[) are also strongly discordant with more detailed analysis presented a few 


years ago (Zamhr et ah, 2010) for the 500 brightest DR5 quasars (Fig. |^. This is 
especially true for so-called population B quasars with FWHM H/? >4000 km s“^ as was 
shown in 2002 (Sulentic et al., 2002| ). Figure 3 in |Marziani et al.| (|2003a|) quantihes the 


dependence of Fell detectivity on S/N and FWHM H/?. Each Shen measure is connected 
to the corresponding Zamhr measure. Taken at face value, automated measures can even 
mislead statistical analyses i.e., it is not that they fail occasionally because of the usual 
problems but they create categories of spurious classes of sources especially for large 
sample sizes. This is shown by the left panel of Fig. where the “outliers” in the optical 
plane of 4DE1 turn out to agree with the main sequence once more accurate measures are 
used. An additional problem is that HellA4686 emission is apparently included as part 
of the Fell blue blend in the Shen et al. ( 2011[ ) data. HellA4686 measurements are not 
reported there. HellA4686 is pernicious leading to overestimation of HellA4686 strength. 


overestimation of Fell template width and claims of an Fell redshift relative to H/3 (Hu 


et al.| 2008 Sulentic et al., 2012). This is shown by Fig. I where the one Fell strong 


source is more straightforwardly interpreted in terms of weak Fell, and strong and very 
broad HellA4686. 

The new analysis might not have been “two steps back” if all quasar spectra were 
the same. Our work from 1989-2001 also emphasized the diversity of the spectroscopic 
(and other) properties of quasars. The concept of two quasar populations (A and B) 
was introduced to further emphasize this point - it is a simplihed analogue of the seven 
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principal spectral types identified in the stellar H-R Diagram. Naturally low S/N spectra 
tend to diminish this diversity leading to the quite natural temptation to indiscriminately 
average many noisy spectra together in order to produce higher S/N composites. But 
without a context with which to average, these composites it confuses rather then clarihes 
the quasar phenomenology and raises an impediment to improved physical models. What 
would a composite of OBAFGKM stellar spectra give us? 

Figure 1 of SH14 allows us to visualize the main El trend but it does not improve 
our understanding of the distribution of sources in the 4DE1 optical plane because the 
Fell blends, present in virtually all quasar spectra, cannot be properly detected, much 
less modeled, in 90-95% of the SDSS spectra. A continuum is likely £t on top of noisy 
Fell emission resulting in a zero, or too low, measure of Fell strength. A censored data 
analysis is appropriate in this context (Sulentic et ah, 2002). 


4. 2002 — Eigenvector Binning and First Explorations of [Olll]A5007 


4.1. Defining spectral types 

We are still speaking about pre-SDSS years. Further exploration/exploitation of 4DE1 
required a larger sample but fewer then 200 quasars with moderate/high S/N spectra 
existed at this early time. While engaged in obtaining new data with 2m class telescopes in 
Spain, Mexico, Italy and Chile we realized that another approach (following OBAFGKM) 
would be to bin the optical plane in order to better contrast quasars along the 4DE1 main 
sequence. Bin sizes were chosen so that all quasars within a given bin were statistically 
indistinguishable. If we were to bin stellar spectra we would bin by spectral type or even 
subtype since indiscriminate averaging will obscure any physical insights. Most of our 
papers showed that tremendous spectral diversity existed in the type-1 quasar population. 
We chose bins of equal size because we have no evidence that the physical drivers of 4DE1 
source occupation correlate with measured properties in a way to justify unequal bins. 
We are not yet in a position to know for sure. |SH14 rediscover binning but with bins of 
unequal size which seems to us a step back. Binning should be predefined rather than 
driven by a distribution in a parameter space dehned using low S/N spectra. Our hrst 


paper showing average spectra of H/? (Sulentic et ah, 2002) revealed a clear change in the 


prohle shape along the sequence with prohles for sources with lowest inferred Eddington 
ratio showing an extra very broad and redshifted component (VBG). This is not to say that 


the binning done by Sulentic et ah (2002) isolated sources that were scattering randomly 


around an average value in each bin: rather, the grid of Marziani et ah (2001) clearly 


shows that within each bin orientation and Eddington ratio trends were expected. 


4-2. Very broad component and very broad line region 

This and other striking differences motivated us however to separate quasars into 
two populations: population A, involving mostly radio-quiet (RQ) high accretors with 
symmetric Lorentz-like H/? prohles, strong optical Fell emission, a ClvA1549 blueshift / 
asymmetry and a soft X-ray excess. Population B sources are lower accretors including 
most of the classical radio-louds quasars (RL). They show composite double-Gaussian 
(BG-I-VBG) H/? prohles, weaker optical Fell emission, unshifted (or redshifted?) ClvA1549 
and no soft X-ray excess. Some quasars show an unusually strong VBG component that 
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dominates the H/3 profile e.g. (PG 1416-129 Sulentic et ah, 2000c). At first the BC + 


VBC decomposition was purely phenomenological: but extreme VBC sources give us a 
clearer insight into the shape of H/3 VBC. In most cases it was possible to reproduce 
H/3 well for the overwhelming majority of population B sources and specihcally the ones 
classihed as AR,R or AR,B (i.e., redward asymmetric, AR, with red or blueshifted peak, 
Sulentic||1989 ). A very broad line region (VBLR) has been postulated since the mid 1980s 
and was initially suggested to be due to optically thin gas in the innermost BLR (Peterson 
and Ferland, 1986 Zheng, 1992; Shields et al.| 1995; Morris and Ward, 1989). Optical 


“thinness” has strong implications for the maximum luminosity associated with a line: if 
the medium is optically thin the intensity of the same recombination line is governed by 
the volume and density of the emitting gas and is not directly related to the luminosity 


of the ionizing continuum (Marziani et ah, 2006). 


Earlier observational dehnitions were based mostly on the ClvA1549 prohle. Broth 


erton et ah (1994a) postulated the existence of a intermediate line region emitting most 
of ClvA1549 of FWHM ^ 2000 km s“^, and a broader component (VBC) to account for 


the extended ClvA1549 wings. Corbin (1995) interpreted the redshift often observed in 


the H/? wings as a gravitational redshift affect the VBC. By 2002-2003 it was clear that 
the prohle analysis of the ClvA1549 prohle was more ambiguous than the one of H/?. 
While H/3 customarily shows a narrower component sharply separated from the broad 
component, the ClvA1549 prohle of Pop. B sources shows a prominent semi-broad core 
that merges smoothly with the line base. We already had suggested an interpretation in 


terms of density stratihed NLR to account for the ClvA1549 core prohle (Marziani et ah 


1996 Sulentic and Marziani, 1999). 


No VBC is apparently present in Fell: we run several tests and we concluded that 
the VBC cannot be as strong in H/3 and, if present, is too weak to be appreciable. The 
BC and VBC decompositions therefore rehects the ionization stratihcation within the 
BLR, simplifying what is probably a continuous trend. Two main regions are identihed: 
a higher ionization VBLR, associated with the VBC (i.e., line base and line wings), and 
a lower one associated with the BC, that emits most of all of Fell. The absence of Fell 
emission can be considered a dehning property of the VBLR. 

It was later shown that the VBC is strong in luminous quasars, and that its luminosity 
cannot be explained by optically thin gas. Presently, the VBLR is understood as inner 
emitting region whose high ionization degree leads to a lower responsivity to continuum 


changes (Korista and Goad, 2004 Goad and Korista, 2014). 


4-3. Blue outliers: linking wind signatures on widely different spatial scales? 

During this year we also began exploration of [Olll]A5007 as a diagnostic of the nar¬ 
row line region. This kind of study requires very high resolution. In fact BG92 expressed 
frustration with attempts to parameterize [Olll] A5007 using standard measures like equiv¬ 
alent width (EW). The reason for this was that even the very good spectra employed in 
BG92 were too low resolution to allow realization that there are two distinct components 
of [Olll] AA4959,5007. We obtained spectra at San Pedro Martir, Calar Alto and KPNO 
and identihed a class of “blue outlier” sources (Zamanov et al., 2002) where the velocity 


separation between the two [Olll]A5007 components (unshifted-l- blueshifted) was largest, 
and the prohle appeared with a blueshift > 250 km s“^ in amplitude. They appear to 
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favor population A quasars (Fig. |^which also show large ClvA1549 bluesliifts, suggesting 
a kinematic linkage between the narrow and broad line regions. Further examples were 
presented the next year ( Marziani et al.| |2003c ), and the existence of large blueshifts and 
semi-broad [Olll]AA4959,5007 emission has been conhrmed in a number of studies (Aoki 

2005 Komossa et al.| 2008; Zhang et ah, 2011). Recent stud- 


et ah, 2005 Bian et ah 


ies even hnd emission extended on galactic scales giving rise to an integrated semibroad 


prohle (Cano-Diaz et al., 2012), although the original results of Zamanov et ah (2002) 


indicated that the semibroad [Olll] AA4959,5007 component originated in a very compact 
NLR, ~ 1 pc in size. At that time and with our data the question mark on the BLR-NLR 
linkage was wise. Nowadays it seems that quasar outflows are linked from the BLR up to 
the circumnuclear regions, where molecular outflows are detected, at least in the extreme 


source Mrk 231 (Feruglio et ah, 2015 Tombesi et ah, 2015), a BAL QSOs that was noted 


in our early analysis because of the abnormally large H/3 FWHM that placed the source 
outside the main sequence of 4DE1 (Marziani et al., 2001 Sulentic et al., 2006). It is 


unfortunate that the binning in SH14 was nonuniform because it damped out some of 
the [Olll]A5007 prohle diversity in their Figure 2. Bin size should not be motivated to 
maximize the number of quasars within a bin but rather to explore the diversity of Type 
1 spectroscopic properties within the parameter space. We have evidence that even our 
already rather hne subdivision should be made even hner. Given the complexity of Fe 
llopt emission and weakness of [Olll] AA4959,5007 in extreme population A sources one 
must be especially cautious interpreting results for these quasars. 

The essential point that emerged in 2002/2003 is that [Olll] AA4959,5007 shows two 
components. The narrower unshifted component shows a strong change along the 4DE1 
main sequence becoming very weak or absent in extreme Pop. A (highest L/LEdd) sources. 
A second semibroad and blueshifted component may be present in all sources but is most 
prominent in extreme population A where it shows blueshifts up to 1000 km s-^ (I Zw 
1 is a famous example of the so-called blue outliers). This is not a novelty of SH14| since 
this trends was discussed in Marziani et al. (2006) on the basis of the previous 2002-2003 
papers. The semibroad component is difficult to study in many sources requiring spectra 
of very high S/N and resolution. In some of these rare spectra [Olll] AA4959,5007 is fully 


resolved into two components (Grupe et al., 1999). A further problem mentioned in some 


of the above references involves slit effects connected with extended emission hne region 
(EELR) contaminating the spectrum. 


5. 2003: an useful, uniform dataset after Boroson and Green (1992) 


5.1. A physical dichotomy between Pop. A and B? 

This was “a very good year” when observations in Mexico, Spain, Italy and Ghile 
(ESO) enabled us to expand our sample to more than 200 low z quasars. This was a big 
deal on the eve of the SDSS. Even better we began to explore in more detail of role of 
black hole mass and Eddington ratio along the quasar main sequence. We were not the 
hrst ( Marziani et al.| 2003c) but the 4DE1 approach allowed us to gain some insight, and 
to suggest that the separation between Pop. A and B occurred as a fairly well dehned 
L/L^dd, ~ 0.1—0.2 for a quasar with log Mbh~ 8. It is possible that this critical Eddington 
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ratio may signal the transition to a slim disk from an optically thick, geometrically thin 
disk ( Abramowicz et ah| 1988). 

Are population A and B simply extreme ends of a main sequence or do they repre¬ 
sent two distinct quasar populations? This remains an open question as far as we are 
concerned, although there is now evidence that the A/B differences could be associated 
with an accretion mode transition, with higher L/L-^dd sources accreting in an advec- 


tion dominated accretion flow (ADAF) (Marziani et ah, 2014, and references therein) 


The gravitational explanation for the large redshift observed at the line base of H/? was 
found to be too large to be consistent with the assumption of predominantly virial motion 
(Marziani et ah, |2003c ). 


5.2. The RL/RQ dichotomy 

It was also the year when more detailed studies of the RQ 


RL dichotomy were ex¬ 


plored in the 4DE1 context hrst with our own atlas 

(Sulentic et al. 

2003 

) and later with 

the 500 brightest SDSS DR5 quasars (iZamhr et al. 

2008 

). The strategy here was to ex- 


plore the distribution along the 4DE1 main sequence of the most unambiguous class of RL 
quasars-those showing classical double-lobe (LD) radio morphology. They clearly show a 
strong concentration at the low LjL-^dd (Population B) end of the main sequence. Almost 
all show a consistent absence of a soft X-ray excess and CIV blueshift. The restricted 
4DE1 parameter space occupation of LD RL sources is perhaps the strongest evidence for 
a physical dichotomy between RL and RQ quasars (Fig. |^ - if 4DE1 parameters reflect 
fundamental aspects of BLR structure and kinematics. The situation is much less clear 
for core-dominated RL - a few are interpreted as preferentially aligned LD sources but 
the majority especially those weaker than the LD sample lower limit (log Pi. 415 ghz ~ 31.6 
ergs s“^IIz) have been attributed to diverse scenarios (e.g. “young” van Breugel et ah 
1984 and “frustrated” Fanti et al. 1995 scenarios). If these scenarios reflect reality they 
imply that some CD sources above (log Li. 415 ghz ~ 31.6 ergs s“^IIz“^) may not be RL in 
the classical sense. They will rise out of the RQ population and many will become rapidly 
quenching flat-spectrum radio sources, or eventually cross the LD boundary producing 
LD morphology (as proposed for the compact steep spectrum (CSS) radio sources). It was 
shown in Zamhr et al. (2008) that CD sources in the range logLi. 4 i 5 GHz ~ 30.0-31.6 ergs 
s“^Hz“^ distribute like RQ quasars and not LD RL. Looking at Fig. 4 and 5 of Zamhr 


et al. (2008) makes obvious that 4DE1 has much to say about the RL-RQ dichotomy. 


6. 2004(-h2007): interpreting the rest-frame UV spectrum and the ClvA1549 
profiles along 4DE1 

Up to this point we were not in a position to explore the other 4DE1 parameters in 
much detail. Myriads of ClvA1549 spectra existed by this time including the very good 
spectra from the Palomar surveys (e.g., Barthel et al., 1990) but we had no way to reliably 


estimate quasar rest frames. Since we knew that ClvA1549 prohle shifts showed large di¬ 
versity along the 4DE1 main sequence we could not proceed. We could not use existing 
studies dBrotherton et ah 1994b) because their methods of ClvA1549 prohle decompo¬ 


sition were incompatible with ours and not physically justihed (Sulentic and Marziani 


1999). We argued that a narrow component of ClvA1549 must be cautiously subtracted 















































and that it was stronger in population B quasars, with semi-broad prohles. Our approach 
was motivated by lack of a well dehned critical density associated with ClvA1549 (instead 
suppressing [Olll]AA4959,5007 at density 10® cm“^. The assumption of a density 
radial trend gave additional arguments that ClvA1549 narrow component (NC) could be 
broader than [Olll] AA4959,5007 by up to 3 times following the simple modelization of 
etzer (1990). ClvA1549 emission it favored at high ionization parameter, with a steep 
decrease toward lower ionization conditions. The emissivity of H/? has instead a much 
flatter dependence on ionization parameter, and its volume integrated prohle is weighted 
toward the outer emitting regions. Hence, the ClvA1549 NC prohle may merge smoothly 
with the BC, while the H/3 NC prohles stands out as an easily separable feature that is dis¬ 
continuous from the BC. Since CivA1549nc is due mainly to the inner NLR, CivA1549nc 
blueshift within a few hundreds km s“^ per second are expected (Zamanov et ah, 2002 


Aoki et al. 

2005; 

Boroson, 

2005 

Bian et al. 

2005 


Komossa et ah, 2008). 


We also gave a recipe for estimating if NC is present and how to subtract it. Obviously 
if one does not distinguish population A and B things look more confusing. If ClvA1549 
NC is present and you do not subtract then you measure ClvA1549 too narrow, too strong 
and probably less shifted. Baskin and Laor (2005) made a NC subtraction although it was 
much lower for sources in common. In population B, where the NC merges smoothly with 
the BC, an NC dehnition is often ambiguous. The most controversial situation indeed 
concerns population B sources. There the BC merges smoothly with the NC prohles, 
making the separation of NC and BC operationally ambiguous. However, not subtracting 
the NC may introduce a signihcant bias lowering Mbh mass estimates, and making them 


fortuitously consistent with the H/3-derived ones (Sulentic et ah, 2007). 


Finally the number of low-z spectra with ClvA1549 spectra in the HST archive passed 
100 enabling us to produce binned spectra ( |Bachev et al. 2004; Sulentic et ah, 2007). Bin¬ 
ning (and rest frame estimation) came from previously obtained matching optical spectra 
of the H/3 region. Now it was very clear that ClvA1549 blueshifts were a Population A 
phenomenon likely associated with a disk wind or outhows in these highest accretors (Fig. 

§■ 

This year also saw the beginning of our study of high redshift quasars in the 4DE1 
context. This required 8m class telescope time because we were forced to follow H/1 -|- 
[Olll]A5007 into the infrared in order to both assign high z sources to the correct bin and 
to obtain a reliable estimate of the rest frame. We used VLT-ISAAC to obtain spectra of 
unprecedented S/N for 52 high 2 ; sources. We could, for the hrst time, observe changes 
in 4DE1 source occupation and we could explore UV surrogates for estimating rest frame 
and black hole mass - lessening the need for IR spectra of large samples. This and much 
more await rediscovery and we certainly hope it will not take another 15 years. 

Equivalent measures for items 1 and 2 (^ have been identihed in the redshifted 


UV spectra of high-redshift quasars (Negrete et al., 2014). In the simplest consideration 


of quasar diversity involving Population A and B sources (described below), virtually 
all multiwavelength measures of quasars show differences (see Table 5 of Sulentic et al. 


2007). 
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7. 2009 


Until the mid 2000s our main results were based on the study of low-z, predominantly 
low luminosity sources. In 2009 we presented the largest instalment of IR spectra of 


Hamburg-ESO intermediate redshift quasars covering the H/3 spectra range (Marziani 


et ah, 2009). In total, considering previous batches (Sulentic et ah, 2004, 2006), we had 


data available for 52 very luminous quasars. Comparison with previous observations at 
lower L allowed as to start the analysis along the L dependent “eigenvector 2.” Perhaps 
the most remarkable result was the detection of the systematic increase in minimum- 
FWHM as a function of luminosity, a trend that is expected if: (1) the H/3 line prohle 
is virially broadened; (2) L/LEdd=l is a physical limit; and (3) the BLR radius scales as 
L oc r" with luminosity, with a ~ 0.5 — 0.65 ( Bentz et al.| 2013; Kaspi et ah, 2005). We 
cannot hnd high luminosity narrow line sources if these assumptions are valid. Negrete 


et ah (2012) and Marziani and Sulentic (2014) have shown the equivalence of a broader- 


lined source with the prototype I Zw 1. This implies that the limit at FWHM = 2000 
km s“^ for dehning Narrow Line Seyfert 1 sources makes sense only at low-luminosity. 
Relaxing this limit would allow to appreciate that there are sources which are analogous 
to local Fell strong NLSyls in terms of physical condition but simply with broader lines. 
A more physically oriented criterion could attempt to isolate sources at L/LEdd—t C>(1), 
which means applying the empirical criterion i?Feii^l (he., to isolate the sources we called 


xA in Marziani and Sulentic 2014). A second intriguing result was the high frequency of 


blueshifted, low equivalent width [Olll]A5007 prohles. Since not all HE sources were high 
L/L-^ddi we suggested that this trend could be associated with a NLR evolutionary effect 
(see also 


8. 2014, or the year of the final rediscovery 

8.1. Results involving [Olu]XX4959,5007 

SH14 establish in a hrmer way the sequence of [Olll]AA4959,5007 prohle behaviour. 


something also described much earlier (Zamanov et al.l 2002 Marziani et ah, 2003b 


2006 Marziani and Sulentic, 2012). Figure supports the trend concerning shifts of 


[Olll]A5007 shown in their Extended Data Figure 2. Its signihcance was discussed in the 
above references and also by Bian et ah ( 2005[ ); Hu et ah (2008); Komossa et ah (2008). 
The panel SFI of SH14 shows the [Olll] AA4959,5007 “Baldwin effect,” i.e., a systematic 


decrease of [Olll]AA4959,5007 EW with L. It has already been stressed (Netzer et ah 


2004 


2006 


trenc 


, and is consistent with the VLT data of Hamburg ESO quasars (Sulentic et ah, 2004 


Marziani et al.l 2009). It is interesting that Shen and Ho (2014) hnd a diheren; 


for the [Olll] AA4959,5007 blue wings - in the sense that wings remains prominent 
also at large luminosity, and shows little or no Baldwin ehect. These hnding may support 
an evolutionary interpretation of El (Dultzin-Hacyan et ah, 2007 Marziani et al.| 2014): 
at the extreme population A end, small Mbh black holes accreting at high rate, with 
evidence of circumnuclear star formation (e.g., Sani et ah , 2010); at the end of population 
B, mostly very massive quasars associated with low accretion rate and old structures like 
extended radio-lobes of Fanaroh-Riley H sources. The evolutionary interpretation seems 
straightforward for RQ quasars. However, while the properties of extreme population A 
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and B are obviously different, and old radio sources are preferentially found at the extreme 


of population B (Zamfir et ah, 2008), it is still not obvious which sources population A 


or B sources may be considered the radio loud progenitors of FRII in the 4DE1 sequence. 


A tantalising possibility is that compact-steep sources (CSS) may evolve into FRII (van 


Breugel et al. 

1984 

Sulentic et al. 

2015) 


8.2. Physical interpretation of the optical El diagrams 

SII14 adopt an approach that is at hrst as surprising as it is unconventional. They 
consider that quasars hosting massive black holes - associated with larger hosts - should 
form in denser environments. They conclude that quasars with higher Rpeii are on average 
less massive because they form in less dense environments. The deduction is intriguing 
and, as presented, has a statistical value: it supports the idea that there is also a trend of 
mass with Rpeii that may lead to a sequence of L/Mbh if the L range is not large. This 
would account for a displacement in the horizontal direction i.e., Rpeii, but not in the 
direction of FWHM. They then conclude that the spread in FWHM is not due to Mbh 
and may largely be due to orientation. As a proof it is incomplete. Their approach also 
follows the hypothesis that it is not possible to have strong Fell in very massive sources. 
However, this is contradicted by our VLT-ISAAC observations of 52 high luminosity ESO- 


Hamburg quasars where we often End strong Fell (Sulentic et ah, 2004, 2006 Marziani 


et ah, 2009) that led to even suspect an “anti-Baldwin effect” for Fe llopt- In our opinion. 


a more convincing proof is offered by at cr* systematically decreasing with increasing Fell 
at hxed quasar luminosity. If Mbh and host galaxy bulge mass are tightly correlated 


(Ferrarese and Merritt 

200C 

Gebhardt et al. 

2000 

), then Fell strength increases with the 

L/M ratio (Sun and Shenf 

2015). 



In order to test the effect of orientation SH14 consider the ratio between black hole 
mass (derived from the Mbh ^ cr* relation) and the virial product / = GMBH,cr/'cFWHM^. 
This is basically the structure factor that should be independent of line width if not 
influenced by viewing angle. A dependence on orientation is however expected from the 
assumption of a flattened, axially symmetric broad line region. They indeed find a strong 
dependence and conclude that orientation strongly affects FWHM H/3. The orientation 
effect has been extensively explored, and a strong dependence conhrmed, using radio-loud 


quasars where orientation can be inferred from radio morphology (Wills and Browne 


1986 


2012 


Sulentic et ah, 2003 Jarvis and McLure 


Zamhr et ah, 2008[ Runnoe et ah 


2013). This rediscovery adds evidence to the hypothesis that orientation is indeec 


vertically displacing sources in the 4DE1 optical plane (as clearly also seen in Marziani 


et ah, 2001). It leaves open the possibility that / is changing due to other effects (for 


example Eddington ratio, although it is reasonable to expect that the effect of L/LEdd is 
not as large as that due to orientation). 

The next conclusion in SH14 is quite surprising: that orientation is driving the change 
in Balmer line prohles between Population B (broader) and Population A (narrower) . A 
signihcant vertical displacement is expected from Mbh (since FWHM oc \/Mbh)- 


manov and Marziani (2002), in a paper that complemented Marziani et ah (2001) showed 


Za- 


the effect of changing Mbh on the optical plane. We also considered the hypothesis that 
orientation was the dominating factor in 4DE1 and that it was driving Pop. A and B 
differences. We concluded that this is very unlikely: there are differences not only in line 
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profiles but also in diagnostic line ratios that cannot be reconciled with an orientation 


effect, and demand a change in physical conditions. The grid in Marziani et ah (2001) - 
derived also considering trends in diagnostic ratios - was computed for a fixed black hole 
mass. It showed that the effect of orientation is most important for sources with weak 
i?Feii- If more luminous sources are included we find a vertical spread associated with 
increasing mass: i.e., large i?Feii and large FWHM. Considering the main effects, a 4DE1 
contextualization must involve at least 3 dimensions in order to distinguish between the 
effects of orientation and mass. And SII14 still results do not provide an estimation of 
viewing angle that is valid for individual sources. 

Accepting at face values the results of Marziani et ah (2001) and SH14| the conclusion 
would be that spectral type Al (i?Feii< 0.5, FWHM < 4000 km s“^) is mainly due to 
face-on population B sources. Indeed an inspection of the sources in spectral bin Al for 
the sample of Marziani et al. (2003c) reveal a composite distribution. Bin Al seems to 
be populated by sources which could be population B oriented pole on (i.e., powerful CD 
radio-loud sources showing a relatively narrow BC and a faint VBC) but also that are the 
apparent extension of A2 with lower i?Feii- The latter type slightly outnumbered the type 


2010), yielding a median Al 


B sources in our samples (including the one of Zamfir et al. 
prohle consistent with a Lorentzian function. 

8.3. 4DE1 parameters 

Lowest EW ClvA1549 sources colored black/blue are favoring the lower end of the 
main sequence in the i?Feii vs FWHM H/? plane. This rediscovered trend in SH14 SFIV 
was previously reported in Bachev et al. (2004) and Sulentic et al. (2007). Note that 


sources in Population bin Al do show relatively large W(ClvA1549). ClvA1549 equiv¬ 
alent width is especially useful for identifying spectral bin type of high redshift quasars 


Negrete et al. (2014) discuss UV identihcation criteria in the 4DE1 context showing that 
50 A is a sufficient criterion for identifying extreme Pop. A sources (bin 
but not Al). ClvA1549 as well as intermediate ionization line 


W(ClvA1549) 
types A2, A3 and A4 
widths and profile shapes distinguish bin Al from Bl. 

Early work resulted in the inclusion of Pgoft as a principal 4DE1 parameter. The pres¬ 
ence of a hard X-ray power-law may represent the best “operational definition” unifying 
all AGN but it was found to show less intrinsic dispersion than the soft X-ray photon 


index (Brandt et al., 1997). There is a clear trend along the main sequence but SH14 can 


only rediscover in their SFVH what we have already learned (Wang et ah, 1996; Sulentic 


et ah, 2000b). It is the availability of X-ray spectra that limits our statistical analyses. 


The main sequence trend is clearest when comparing mean/median measures for the up¬ 
per (population A) and lower (population B) ends. Spectral types A2-A4 show the largest 
soft X-ray excesses and Pop. B only show a hard power/law. A comparison of all available 
(up to 2006) XMM spectra (n = 20 - 40) for Pop. A and B sources ([Sulentic et al.[ |2008[) 


confirms the differences originally reported in Sulentic et al. (2000b) in the sense that only 


Pop. A quasars show a soft excess. Note that the papers providing XMM data for this 
study consistently report no trends. SH14 show more Chandra/XMM data points than 
were available in 2006 but it is unclear what filters have been applied to existing X-ray 
data (e.g number of hard/soft X-ray photons detected). 
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8.4- RL/RQ dichotomy 

There has been considerable work seeking connections between 4D Eigenvector 1 pa¬ 
rameters and radio properties both to nse radio-loud sources as orientation indicators and 


to explore the RQ-RL dichotomy 

(Wills and Browne 

1986 

|Marziani et al. 

2001 

Rokaki 

et al. 

2003 

Sulentic et al. 

2003; 

Zamhr et al. 

2008 

)• 

SH14 

wisely distinguish between 


core and lobe dominant RL sources in their SFIX. The problem with this “rediscovery” 
and a step backwards comes from their adopted definition of radio-loudness. A widely 


used criterion sees as RL all sources with radio/optical flux ratio i?K>10 (Kellermann 
et ah, 1989). Our own work (Sulentic et ah, 2003| Zamfir et ah, 2008) has focussed on 
lobe-dominated (LD) RL sources as the parent population of classical RL quasars. LD 
sources show a well defined lower limit in radio power (i?K>70) and a restricted occu¬ 
pation along the upper end of the 4DE1 main sequence. CD sources do not show these 
properties and those between i?K=10-70 distribute along the main sequence like the RQ 
quasar majority. Inclusion of these sources in a RL sample will tend to obscure differences 
(and dichotomy) between RQ and RL quasars. Perhaps ironically, a larger sample of the 
brightest 150 LD quasars from SDSS DR9 conhrm our original lower radio power limit 
for these sources as well as their restricted 4DE1 population B occupation. 


9. Important open issues 

The previous sections outline an undeniable progress in the understanding of quasar 
multifrequency properties, with a successful contextualisation in eigenvector 1 based 
schemes. We think that there are however some open issues that should be faced openly 
to set eigenvector 1 studies on a hrmer ground. The hrst one is undeniably the forma¬ 
tion of Fe llopt and FeiIuv lines. Current photoionization calculations do not allow to 
reproduce sources with Rpeii^ 1, and it is still unclear whether this is a shortcoming of 
photoionization codes (that being based on a mean escape probability formalism for the 
treatment of radiation transfer, are not suited to consider low ionization line formation in 
the extended partially ionised zone where absorption processes are mainly non-local). On 
the one hand a study of other low ionization lines like the Call IR triplet indicate that a 
dense, low ionization photoionizing medium can reproduce their total emission even in the 


most luminous sources (Matsuoka et ah, 2007 Loli Martinez-Aldama et ah, 2015). On the 


other hand, extreme Pop. A sources with Rpeii^ 1 are still poorly understood. They are 
most likely affected by a strong wind that may lead to shielding of the continuum (e.g.. 


Elvis 

2000 

fLeighly and Moore 

2004 

Leighly et al. 

2007 

) or to mechanical heating of the 

line emitting gas. The relatively old paper by 

Marziani et al. 

(2001 

) estimated a decrease 


in U with increasing L/L-^dd, but connect U and Fell prominence on an empirical relation 
and on crude scaling assumptions. Another aspect, that is severely constrained by data 
availability, is the interpretation of the Tgoft that may complicated by the presence of the 


warm absorber (Done et ah, 2012 Chakravorty et ah, 2012). Last, we have to mention 


that not all workers agree that the ClvA1549 blueshift can be interpreted in terms of an 


outflow (Gaskell and Goosmann, 2008), and that the best approach to extract a virial 


broadening indicator from the line width of H/3, especially in Pop. B, is still a debated 
issue (e.g., Gollin et ah, 2006 Marziani and Sulentic, 2012 Shen, 2013). In this context. 
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it remains important to ascertain the nature of the redward displacement of the VBC in 
Pop. B sources (Sulentic et ah, in preparation). 


10. Conclusion 


SH14 conhrm several empirical trends that were discovered earlier. Among them the 
ones involving [Olll] AA4959,5007, distributional differences of RL/RQ sources in 4DE1, 
H/3 prohle properties as well as soft X-ray and CIV measures. SH14 adds an empirical 
verihcation that orientation effects matter also for RQ sources - something missing from 
earlier papers that used the FWHM dependence on orientation derived from RL samples. 
It does not clarify the physical factors driving the location of sources in the 4DE1 optical 
plane. Only a multidimensional approach involving at least Mbh, L/L^dd and orientation 
(cf. Laor, |2000 ; [Marziani et ah 2001 Zamanov et al.| 2002) can lead to a hnal clarihcation 
of these issues and successful development of a physical model. 
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Figure 1: Simplified version of Fig. 4 in Marziani et al. (2001). The 4DE1 optical plane shows sources 
from the Sulentic et al. (2000b) sample with a superposed grid of theoretical values of viewing angle f (10 
<i< 40) and L/M, expressed in solar values, for 3.1 < logM < 4.5 (implying log T/TEdd~ 0), at steps 
of Slog L/M = 0.1. A value of logM « 8 in solar units was assumed for U. Filled circles and squares 


represent RQ and RL sources respectively in the Sulentic et al. 

(2000b 

) sample. Boxes identify BAL 

QSOs (2 of them, Mark 231 and IRAS 0759+651 not in 

Sulentic et al. 

2000b 

). The red spot mark the 

average for the strong Fell emitters of 

Lipari et al. 

(1993 

. Error bars in 

the upper right corner indicate 
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Figure 2: Left: a comparison between Shen et al. (2011) and Zamfir et al. (2010) in the optical plane of 
4DE1 involving the width of H/3 and the strength of optical Fell. The blue lines connect [Zamfir et al.| 
(2010) (green stars) with Shen et al. (2011) measurements for the same objects. Right: A2 and part of 
adjacent bins in the 4DE1 optical plane 
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Figure 3: Top: continuum placement for source SDSS J034106.75+004609.9. The black line traces the 
featureless continuum; the green line traces the total continuum that includes a contribution from the 
host galaxy (whose presence is made manifest by the Mg Ib band visible in the spectrum. There is no 
significant Fell emission. The magenta line traces the cleaned H/3 broad profile. Bottom: Fl/3 profile 
analysis. The excess on the blue side of H/3 can be explained by significant HeiiA4686 VBC emission. H/3 
and Hy are consistently fit with the same BC and VBC parameters (the BC and VBC (red) decomposition 
is shown for H/3 only). 
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Figure 4: The 4DE1 optical plane involving the width of H/3 and the strength of optical Fell. Horizontal 
line marks the boundary between population A (lower) and B (upper) sources. Source occupation is 


shown for the 470 brightest SDSS-DR5 quasars (z<0.75) with highest s/n SDSS spectra Zamfir et al. 


(2008). Grey and red dots represent RQ and lobe dominated LD RL quasars. The latter show a strong 
preference for the Population B zone. 
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Figure 5: Same optical plane as the previous Figure, for all low z quasars with measurable HST/FOS UV 
spectra (Sulentic et al. 20071. Different symbols represent the amplitude of the CIV 1549 blueshift at 


half maximum c(^). Large blue filled circles involve sources with largest CIV blueshifts which strongly 
favor the Population A zone. Large open red circles represent sources with largest CivA1549 redshift and 


grey squares those with no significant line shift (Sulentic et al. 20071. 
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Figure 6: Average shift of [Oiii] A5007 (negative is blueshift) for spectral types in the optical plane of El. 
Numbers in parentheses indicate sources in each bin, with the dataset available to us in late 2001. 
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